Effects of Random Density Fluctuations on 
Matter- Enhanced Neutrino Flavor Transitions 
in Supernovae and Implications for Supernova 
^5 ; Dynamics and Nucleosynthesis 

m 

<N : F. N. Loreti, Y.-Z. Qian, 

G. M. Fuller? and A. B. Balantekint 
\q ; Institute for Nuclear Theory, 

University of Washington, Box 351550, Seattle, WA 98195 

oo 
O 
in 

G\ ■ Abstract 



Oh 



We calculate the effects of random density fluctuations on two- 
neutrino flavor transformations [y r (ji) ^ v e) i n the post-core-bounce 
5^ ■ supernova environment. In particular, we follow numerically the fla- 

vor evolution of neutrino states propagating through a stochastic field 
of density fluctuations. We examine the approach to neutrino flavor 
depolarization, and study the effects of this phenomenon in both the 
^ . early shock reheating epoch and the later r-process nucleosynthesis 

epoch. Our results suggest that significant fluctuation-induced neu- 
trino flavor depolarization effects occur in these environments only 
when the zero-order (without density fluctautions) evolution of the 
neutrino states includes adiabatic propagation through resonances 
(mass level crossings). In the shock reheating epoch, depolarization 
effects from fluctuations with amplitudes larger than 0.05% of the lo- 
cal matter density can cause an increase in the heating rate of the 
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material behind the shock compared to the case with no neutrino fla- 
vor transformation - but this corresponds to a significant decrease in 
this quantity relative to the case with adiabatic neutrino transforma- 
tion. If r-process nucleosynthesis is to occur during the late stages 
of supernova evolution, then the requirement of neutron-rich condi- 
tions excludes a region of neutrino mass-squared difference and vac- 
uum mixing angle (5m 2 , sin 2 20) parameter space for neutrino flavor 
transformation. We find that in the presence of stochastic fluctua- 
tions, this excluded region is not significantly altered even for random 
fluctuations with an amplitude of 1% of the local matter density. 
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1 INTRODUCTION 



In this paper we investigate matter-enhanced neutrino flavor transformations, 
the MSW (Mikheyev-Smirnov-Wolfenstein) effect [1,2], for the case of a mat- 
ter density consisting of a smooth component and a fluctuating part modeled 
by Gaussian colored noise. We will consider the effects of fluctuation-induced 
neutrino flavor decoherence during two post-core-bounce epochs of supernova 
evolution: (1) the shock reheating epoch at TPB (time post bounce) < 1 s; 
and (2) the hot bubble r-process nucleosynthesis epoch at TPB rs 3-16 s. 

The neutrino heating of the stalled shock has been examined by Fuller 
et al. [3] and Wilson and Mayle [4] including neutrino flavor-mixing effects. 
Flavor mixing between a light v e and or u T with a cosmologically signif- 
icant mass in the range of 10-100 eV can increase the supernova explosion 
energy by up to 60%. The increase in the explosion energy is a result of 
the greater temperature of and v T compared to that of v e . Because of 
their charged current interactions with nucleons, v e remain in thermal equi- 
librium with the matter to lower densities and temperatures than and v T . 
Thus, although the neutrinos have approximately the same luminosities, fla- 
vor transformation between or v T and v e on their way to the stalled shock 
after thermal decoupling will increase the average energy of v e , resulting in 
more heat liberated behind the shock than the case without neutrino flavor 
transformation. 

This increase in the average energy of v e by the MSW transformation also 
affects the possible r-process nucleosynthesis of heavy elements in supernovae 
[5]. In the absence of flavor transitions, v e have a higher temperature than 
v e . This is due to the fact that the neutron-rich surface of the proto neutron 
star presents a greater opacity to v e than u e . Therefore, v e decouple deeper 
in the proto neutron star than u e , and hence, are hotter. Due to their lack of 
charged current interactions (at these energies) with nucleons, and v T are 
even hotter than V e . As pointed out by Qian et al. [5], an MSW transition 
can result in a proton production rate which is greater than the neutron 
production rate at the proposed r-process site. By requiring neutron-rich 
conditions at the r-process site, Qian et al. were able to put limits on the 
vacuum mass-squared difference and the vacuum mixing angle between u e 
and or v T . Since and v T have identical energy spectra, we will hereafter 
consider only MSW transitions between v e and v T for convenience, although 
all our results could equally well apply to MSW transitions between v e and 
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A general, semi- quantitative approach to neutrino oscillations in inho- 
mogeneous matter was developed in Ref. [6]. Studies of matter density 
fluctuations which are not random, but harmonic [7,8], or occur as a jump- 
like change in the solar density [9], are available in the literature. Inhomo- 
geneities in the velocity field of material can mimic the effects of density 
inhomogeneities on neutrino flavor oscillations [8]. Velocity inhomogeneity 
effects become important if the characteristic velocity is near the speed of 
light [8]. Noisy mixing of matter could also mimic a fluctuating matter den- 
sity. 

A priori, fluctuations in the matter density may be well approximated 
by random noise (which averages to zero) added to the average value of 
the density. In Ref. [10], a differential equation for the averaged survival 
probability was derived for the case in which the random noise was taken 
to be a delta-correlated Gaussian distribution. It was also shown that if the 
correlation length of the matter density fluctuations was small compared to 
the neutrino oscillation length at resonance, one obtained the same result 
as for the case of a delta-correlated Gaussian. Here, we consider the more 
realistic case of colored noise [11]. We find that the random fluctuations 
have the largest effect on the flavor transition when the correlation length 
is on the order of the neutrino oscillation length at resonance, and that the 
fluctuations, on average, have the effect of suppressing the flavor transition. 
This result is in agreement with more qualitative arguments regarding flavor 
transition in an inhomogeneous distribution of matter [6]. 

Suppressing the degree of flavor transition has the effect of lowering the 
average energy of v e compared with the average energy they would have 
had for an MSW transition in the absence of the noise. We will show that 
neutrino flavor decoherence caused by random fluctuations has little effect 
on the part of the neutrino mixing parameter space (5m 2 , sin 2 29) previously 
excluded by considerations of MSW neutrino flavor transformation in the 
hot bubble r-process region of the supernova environment. However, we find 
that neutrino flavor decoherence can have significant effects on the neutrino 
heating rate and the electron fraction during the shock reheating epoch. 

In Sec. II, we develop the equations for two-neutrino flavor conversion 
in the presence of multiplicative colored noise. We present a differential 
equation for the averaged transition probability, which is valid when the 
correlation length is small compared to the width of the resonance region. 
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Section III presents our results for the effects of random fluctuations on 
neutrino heating of the shock and r-process nucleosynthesis in the neutrino- 
heated supernova ejecta. In Sec. IV, we discuss our results and present our 
conclusions. 

2 NEUTRINO FLAVOR EVOLUTION IN THE 
PRESENCE OF COLORED NOISE 

In this section, we discuss how an initially pure neutrino flavor state prop- 
agating through a stochastic field of density fluctuations can evolve into a 
mixed ensemble of neutrino flavors. We call this process neutrino flavor de- 
coherence (or flavor depolarization). Here we consider two-neutrino mixing 
and we assume that there is the usual unitary transformation between flavor 
eigenstates (e.g., \v e ) and \v T )) and mass eigenstates (\v>i) and l^))'- 

\u e ) = cos0|i/i) + sin 0| 1*2), (1) 

\u T ) = — sin 0| 1/1) + cos6\u2), (2) 

where 9 is the vacuum mixing angle. Complete flavor decoherence would 
occur if a neutrino emitted from the neutrino sphere in an initially pure 
flavor state, for example \v e ), becomes a mixed state of 50% \v e ) and 50% 
\v T ) (both at the original neutrino energy) after propagating through a region 
where density fluctuations exist. 

Our goal in this section is to find which fluctuation characteristics (e.g., 
root-mean-square amplitude, power spectrum, etc.) would be required to 
cause significant deviations in neutrino flavor evolution from that predicted 
by conventional MSW studies with a smooth density run in the post-core- 
bounce supernova environment. We begin with a discussion of the general 
flavor evolution problem and the way in which fluctuations can be charac- 
terized. 

The density matrix of two-neutrino flavor evolution obeys the equation 

i±P=[H,p], (3) 

where 

^=($r))® ( ° :(r) ' <(r)) ' (4) 
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with a e and a T the probability amplitudes for the neutrino to be u e and u T , 
respectively, and where the Hamiltonian is given by |l J 

/ Sm 2 1 \ /Sm 2 \ 

H= cos 26 + -j=G F (N e (r) + iV;(r))J a 2 + [^^29 j a x . (5) 

In Eq. (5), 5m 2 is the vacuum neutrino mass-squared difference, E is the 
neutrino energy, N e and iVJ are the averaged and randomly fluctuating parts 
of the electron number density, respectively, and a x and a z are the Pauli 
matrices. 

For colored noise, we take the ensemble averages of the randomly fluctu- 
ating part of the electron density to be given by 

(K(r)) = 0, (6) 
(N r e (r)N:(r')) = (3 2 N e (r) N e (r') exp(-|r - r'\/r c ), (7) 

with the averages of all odd products vanishing and all higher even prod- 
ucts given by all possible independent products of two-body correlations 
(i.e., the fluctuations are Gaussian). For example, if we define fu... = 
(A^J(ri)iVg (r 2 ) • • •), the average of a product of four would be /1234 = /12/34 + 
/13/24 + fuf23- In Eq. (7), f3 is the ratio of the root-mean-square fluctuation 
to the local density, and r c is the correlation length. 

Clearly, the detailed evolution of the position (time) dependent ampli- 
tudes a e (r) and a T (r) through regions which include MSW resonances (mass 
level crossings) is quite complicated in the presence of fluctuations and so 
necessitates numerical treatment. However, before we proceed to the de- 
scription of the numerical calculation, it is advantageous to define a few 
quantities which we shall later employ. First, we introduce the resonance 
width as 

5r = 2Ht&n26, (8) 

where the effective weak charge density scale height is H = \d In N e (r) / dr\~ l 
[1]. We define A = 5m 2 /2E, where Sm 2 is the difference of the squares of the 
vacuum neutrino mass eigenvalues. With this notation, resonance (neutrino 
mass level crossing) occurs at the position where A cos 2$ = ^/2GfN £ is 
satisfied. 

We now proceed to a detailed description of how we follow numerically 
the time evolution of the neutrino density matrix [i.e., Eq. (3)]. Provided one 
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has access to a random number generator which generates Gaussian deviates, 
Eq. (3) can be integrated by the method of Ref. [11] with the probabilities 
calculated and averaged. However, such a process can be quite time con- 
suming. For the case in which the correlation length is small compared with 
the width of the resonance region, a "ladder" approximation can be made 
allowing one to obtain a differential equation for the averaged probabilities. 

In order to obtain this approximation, one can transform p to the inter- 
action picture, express pi as an iterative expansion and explicitly perform 
the averaging to obtain [10] 

1 00 r 

( Pl (r)) = - (1 + a,) + ^(-1)" / dR(2n)F(2n)C(2n), (9) 

A n=l J 



where 



/rr rri rri n -\ 

dR(2n) = dn dr 2 --- dr 2n , (10) 

JO J o Jo 

F{2n) = /l23-(2n-l)(2n), (H) 

C(2n) = [M(n), [M(r 2 ), [• • • , [M(r 2n ), a z \ • • •], (12) 

and M = (G F /2y/2)Ufa z U , with U Q being the propagator of the flavor 
evolution in the absence of fluctuations. The averaged density matrix is then 
given by (p) = Uq{pi)Uq. From the physical fact that the fluctuations should 
not affect the evolution far from resonance, one can restrict the integration 
limits in Eq. (10) to the resonance region. Since the matrix elements of 
C(2n) are always < 1, we evaluate the integrals of the terms in F(2n) and 
compare the relative size of these terms. When r c <C 5r, the leading term 
in the integral of F(2n) is of order (r c 5r) n . This contribution comes from 
the term in which the subscripts are in the order of the nested integrals [i.e., 
(/12/34 • • ' /(2n-i)(2n))]- Terms which are out of this order by s interchanges are 
of order (5r) n ~ s T™ +s when r c <^ Sr. Note that because of the absolute value 
signs in Eq. (7), f, tj = f jh so that f 13 f 2 4 • • • /(2n-i)(2n) is out of order by one 
interchange and /13/45/26 ' ' ' /(2n-i)(2n) is °ut of order by two interchanges. 
Therefore, the integral of all terms in F(2n) not in the order of the nested 
integrals will be of order r c /5r or smaller than the integral of the term in 
F(2n) which is in the order of the nested integrals when r c <^ 5r. Therefore, 
we approximate Eq. (9) by keeping only the largest term in each F(2n). One 
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then obtains an infinite series which is the iterative expansion of 

(pi(r)) = hi + a z )- fdn f dr 2 (5 2 N e {r 1 )N e {r 2 )eM-{ri-r2)lr c ] 
2 Jo Jo 

x [M^OJM^),^^))]], (13) 

where we have returned to the full integration limits since the contribution 
from positions outside the resonance region is small. Equation (13) could also 
have been obtained by making the assumption (N^N^'pi) ~ (N^N^')^). 
Since the maximal reduction in the transition probability occurs for cases in 
which the correlation length is approximately the neutrino oscillation length 
at resonance divided by n (i.e., r c ~ L res /7r = 4E/5m 2 sin 28), the above ap- 
proximation should be valid for adiabatic transitions for which the oscillation 
length at resonance is much smaller than the width of the resonance region. 

As an example of the approximation, we calculate and present the sur- 
vival probability as a function of the correlation length for a v T traveling 
through the supernova at TPB ps 3 s in Fig. 1(a). The rms value of the noise 
is taken to be 1% of the local electron number density, and we choose the 
following parameters: 5m 2 = 10 2 eV 2 , sin 2 29 = 1(T 3 , and E = 33 MeV. The 
correlation length of the random fluctuations varies from r c = (0.25-24)r°, 
where = L xcs /t^ = (8m 2 sin 26/4E)~ 1 = 832.6 cm. For these parame- 
ters, the neutrino flavor evolution in the absence of noise is highly adiabatic 
(L rcs /5r ~ 0.03). In Fig. 1, the smooth curve is the approximation and the 
jagged line is the simulation of Eq. (9) utilizing numerical averaging. The 
survival probability in the absence of fluctuations is not shown in Fig. 1(a), 
being ~ 0. The agreement between the approximation and the simulation is 
quite good. We also plot in Fig. 1(b) the v T survival probability as a func- 
tion of the correlation length for the same parameters as in Fig. 1(a), except 
with sin 2 26 l = 10~ 5 and the corresponding r c ° = 8326.0 cm. In Fig. 1(b), 
the horizontal line is the survival probability in the absence of fluctuations. 
One observes that at small values of the correlation length [r c ~ (0-5)r°], the 
survival probability from the simulation is larger than that calculated from 
the approximation, although they are not greatly different. 

The reason the approximation works as well as it does is that for both 
t c ^> L rcs and r c <C L res , the approximation and the numerically-averaged 
simulation approach the probability in the absence of fluctuations. As dis- 
cussed in Ref. [10], the parameter 7 = l/2G 2 F {(N^) 2 )T c Sr governs the size of 
the effect of the fluctuations when r c <^ L res . When 7 is large, the transition 
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probability is heavily suppressed, where as for 7 <C 1, the fluctuations have 
little effect. Therefore, a small value for r c will result in only a small change 
in the probability. As r c goes to infinity, 

2n 

lim o F(2n) = (2n - l)!!,3 2n [] ^e(n) 

= -7™/ rfxexp[-x 2 /(2/5 2 )]x 2 ™n^(r,), (14) 

and Eq. (9) can be summed and transformed back from the interaction 
picture to give 

Jim (p(r)> = -j=p J^dxexp[-x 2 /(2^)}p(r,x), (15) 

where p(r, x) is the density matrix calculated using Eq. (3) for the electron 
density (1 + x)N e (r). Equation (15) provides a simple physical picture for 
the averaged density matrix in the limit of very large correlation lengths: one 
simply calculates the density matrix for the electron density with a "frozen" 
fluctuation, (1 + x)N e (r), then averages over all such fluctuations with a 
Gaussian weight. Large correlation lengths imply that fluctuations at many 
different locations are coupled. Indeed, an averaging such as given in Eq. 
(15) is common in other physical situations when many coupled channels are 
present. For example, in multidimensional dissipative quantum tunnelling, 
barrier transmission probability is given by a similar formula when the num- 
ber of channels gets very large [12]. Note that to first order in x, the change 
N e (r) — > (1 + x)N e (r) in the functional form of the density will affect the 
survival probability via a change in the slope of the density at resonance. For 
slowly changing slopes and small f3, the change in the survival probability 
should be proportional to j3 2 , and therefore be quite small. 

3 EFFECTS ON THE r-PROCESS AND 
SHOCK REHEATING 

As demonstrated by Fig. 1(b), even a 1% fluctuation in the matter density 
(a value which is probably unrealistically large for the rather quiescent TPB 
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> 3 s epoch) does little in reducing the transition probability (or increasing 
the survival probability) for non-adiabatic evolution. As shown in Ref. [5], 
the neutrino-heated supernova ejecta must have a neutron excess (Y e < 0.5) 
in order for any r-process nucleosynthesis to be produced. In Ref. [5], it was 
shown that the electron fraction Y e is approximately given by 

Ye w i + KjKn w i + (E P ,)/(Eu.y (16) 

where A PeP and X VeTl are the reaction rates for 

v e +P — >• e + + n, (17) 
v e + n — > e~+p, (18) 

respectively, and (E Ve ) and (-Ep e ) are the average energy for u e and P e , re- 
spectively. Since at this epoch the average neutrino energies are (E Ve ) ~ 11 
MeV, (£p e ) « 16 MeV, and (E Vt ) fx 25 MeV, a substantial conversion of 
z/ T into z/ e can decrease the ratio (E Pe )/(E Ue ) below one and thereby make 
conditions impossible for the r-process. The average v e energy after an MSW 
transition is approximately P{y e — > u e ){E Ue ) + P{y T — > u e ){E Ui ). To obtain 
neutron-rich conditions, one must have P{y e — > u e ) > 64%. It only requires 
about 30% efficiency in flavor conversion to drive the neutrino-heated super- 
nova ejecta too proton rich for r-process nucleosynthesis. Consider simple 
estimates for Y e using Eq. (16) for four cases: 

(1) No flavor mixing. In this case, we have Y e as 1/(1 + (E Pe ) / (E Ue )) ps 
1/(1 + 16/11) ai 0.41, so the material is neutron rich, in good agreement with 
detailed supernova model calculations [5]. 

(2) Full flavor conversion. In this case, we have (E Ve ) # (E Ut ), so that 
Y e Pd 1/(1 + (E Pe )/(E UT )) Pd 1/(1 + 16/25) ^ 0.61, and the material is very 
proton rich. 

(3) Complete neutrino flavor depolarization. In this case, there is 50% 
flavor conversion, so (E Ue ) — > {(E Ue ) + (E Ut ))/2 fx 18 MeV, which implies 
that Y e ^ 1/(1 + (E Pe ) I (E' Ue )) Pi 0.53, too proton rich for r-process nucle- 
osynthesis. 

In fact, we must have Y e < 0.5 to get any r-process nucleosynthesis. 
This is a conservative limit, since a good r-process requires Y e < 0.45. Now 
consider a fourth case, 
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(4) Partial flavor depolarization. For 35% conversion of v T into v e , we have 
Y e ps 0.5. For 30% conversion of v T into u e , we have Y e ps 0.49, realistically too 
large to give acceptable r-process nucleosynthesis. In fact, to get Y e < 0.45 
we need to demand that there had been less than 15% flavor conversion. 

Consider some condition along the Y e = 0.5 line in Fig. 2 of Ref. [5]. For 
example, consider E u = 25 MeV and 5m 2 ps 900 eV 2 , for which the density 
scale height at resonance is H ps 0.5 km. This value of 5m 2 corresponds to 
sin 2 29 PS 4 x 10~ 6 on the Y e = 0.5 line. In this case, the conversion probabil- 
ity is P(v T -> u e ) « 1 - exp{-0.04 (5m 2 /eV 2 ) (MeV / E u ) (if /cm) sin 2 29} ^ 
25%. The neutrino energies around 25 MeV are the most important in terms 
of leverage on Y e . Note that fluctuation-induced depolarization at a level of 
50% could conceivably produce more conversion than MSW transformation 
at E u = 25 MeV over a considerable part of the region to the right of the 
Y e = 0.5 line in Fig. 2 of Ref. [5]. As outlined above, greater than 30% flavor 
conversion will always drive the material too proton rich for r-process. 

One can conclude that the random fluctuations will have a quite minor 
effect on the neutrino mixing parameters constrained by r-process nucle- 
osynthesis because the absolute increase of the survival probability rapidly 
diminishes with increasing survival probability. To illustrate this, we plot 
the survival probability as a function of energy for v T with fluctuations of 
1% and 0.5% of the local matter density in Figs. 2(a) and 2(b). We choose 
the parameters of Fig. 1(b) so that the evolution is nonadiabatic, and the 
neutrino energy is chosen to maximize the differential capture rate [see Eq. 
(23)]. In Fig. 2, the jagged lines are the solution using Eq. (9), the solid line 
is the approximate solution, and the dashed line is the survival probability 
in the absence of fluctuations. One observes that the approximate solution 
does fairly well in reproducing the simulation. There is a small increase of 
about 0.08 in the survival probability in Fig. 2(a) where the fluctuations are 
1% of the local density, but only a very small increase is obtained for 0.5% 
noise in Fig. 2(b). If one calculates Y e for these parameters in the absence 
of fluctuations, one obtains Y e = 0.5 indicating that this point lies on the 
boundary of the excluded region of the MSW parameter space [5]. The in- 
clusion of noise at the 1% level would decrease Y e by about 1%. Random 
fluctuations of 0.5% as in Fig. 2(b) would give a decrease of about 0.4%. For 
noise at the 0.1% level, there will essentially be no change in the excluded 
region. 
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Note that if the neutrino flavor transition is adiabatic, the fluctuations 
can increase the survival probability from approximately zero to one half. 
The shock reheating epoch occurs at approximately TPB ~ 0.15 s, and the 
relevant scale heights are larger than those at TPB > 3 s, which implies a 
larger resonance width for comparable values of 5m 2 and sin 2 29. Hence, 7 
may remain ~ 1 while ((iVJ) 2 ) is reduced allowing one to obtain a sizable 
effect for adiabatic transitions as shown in Fig. 1(a) for TPB k 3 s with a 
smaller rms value of the noise. In Fig. 3 we show the survival probability as 
a function of neutrino energy for the values 5m 2 = 10 3 eV 2 , sin 2 29 = 10~ 6 , 
r c = 2.762 x 10 3 cm. The rms fluctuations are 0.05% and 0.1% of the local 
density in Figs. 3(a) and 3(b), respectively. In both cases there is a large 
region where the survival probability is increased from the value of essentially 
zero in the absence of fluctuations. The approximate solution seems to be 
better in Fig. 3(b) where the rms value of the fluctuations is larger, but is 
reasonably close to the numerical simulation in both cases. For the case of 
an adiabatic transition between the more energetic v T and less energetic z/ e , 
the heating rate can be increased by (30-60)% [3]. Therefore, the reduction 
of the transition probability by random fluctuation effects for an a priori 
reasonable rms value of the fluctuations [(0.05-0.1)%] could be important. 

We can estimate the decrease in the heating rate due to the fluctuations 
over that from the case of an adiabatic MSW transition as follows. In the 
absence of neutrino flavor conversion, the heating rate per nucleon is approx- 
imately given by [13] 

L v f™Elc vN dE v /{exp(Eu/T„) + l) _ [°° d^{E v ,T v ) 
e ^4yrr 2 / °° ^^/(exp^/T,) + 1) " Jo ^ uN dE v 

(19) 

where 

a^9.6xl0- 44 (^) 2 cm 2 , (20) 

is the absorption cross section, L v is the neutrino luminosity (assumed equal 
for all species), T v is the neutrino temperature, d(p is the differential neutrino 
flux with respect to neutrino energy, and E e is the energy of the produced 
electron or positron in Eqs. (17) and (18). 

The total heating rate combining contributions from both v e and v e is 

^tot = Yn^Ueti ~l" ^pCp e p, (21) 
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where Y n and Y p are the number fractions of free neutrons and protons, 
respectively, with Y n + Y p « 1. Of course, Y n and Y p will be set by the same 
weak interactions responsible for heating the shock. In the region where 
neutrino heating dominates, we can assume 

£ « T*, (22) 

where \ u n is the neutrino capture rate on free nucleons and is given by 

Kn ~ (~ * V N d4 ^ v i : Tv) dE v . (23) 

An MSW transition between z/ e and u T will increase the heating rate by 
a factor of 

•/ Y'e' + F'e- 

fc tot _ 1 n c u e n ^ 1 p c i' e p / 9 .n 
e ~ Ye 4-Yf ' ^ ' 

where e' UeH is given by 

+ Hi - P(K,T e ,P)]E v v vIf Wf' Tvr) dE v . (25) 
We take and K' to be given by Eq. (22) with 



K e n = r P(E U ,T C ,P)<T UN - 
J 



d(p(E u ,T Ue 



-dE v 



/o dhj v 
+ Hi - PiE^M^uN ^/f^ dE v . (26) 

In the above equations, P(E u ,t c , f3) is the survival probability of z/ T (or z/ e ). 

We have calculated Eq. (24) for the cases of random fluctuations with 
an rms value of 0.05% and 0.1% of the local matter density, and for the case 
without fluctuations. We present the results in Table I. For each set of values 
for (5m 2 , sin 2 26 l ), we have taken r c = (5m 2 sin 2#/4i?p eak )~ 1 , where E peak = 
5T„ T , the energy for which the integrand is approximately maximized. We 
have taken the temperatures to be T Ve = T Ve 5 MeV and T„ T 7 MeV. The 
values in Table I were calculated using the approximation in Eq. (13) rather 
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than the "exact" numerical simulation, since we are interested in estimating 
the size of the difference in neutrino flavor conversion efficiency for the cases 
of smooth and noisy density distributions, not on a particular model of the 
fluctuations which may or may not obey Eqs. (6) and (7). From Table I, 
one sees that the noise can reduce the heating rate by up to 45% from the 
adiabatic MSW case. 

4 DISCUSSION AND CONCLUSIONS 

To our knowledge, there is no consensus on the size of possible matter den- 
sity fluctuations in post-core-bounce supernovae. Although we have taken 
the rms size of the fluctuations to be a constant fraction of the local matter 
density, it may very well increase or decrease with decreasing density, or the 
size could depend on the distance from the shock. Similarly, the correla- 
tion lengths we have used were chosen to give the maximal reduction in the 
MSW transition probability and were independent of the density and the dis- 
tance from the shock. If, for instance, the correlation length increased with 
decreasing density, the effect of the noise in reducing the MSW transition 
probability could be enhanced since for fixed Sm 2 and sin 2 29 the correlation 
length of maximal effect (r c ~ L rcs ) varies linearly with the neutrino energy. 
Our intention in this paper is to establish the maximal effect density fluctu- 
ations could have if they can be well approximated by random noise added 
to the average density. 

For an rms fluctuation value of 1%, the addition of the noise will have 
a slight effect on the r-process nucleosynthesis as compared with the MSW 
effect in the absence of noise. This can be traced to the fact that demanding 
neutron-rich conditions at the site of the r-process eliminates all but the 
nonadiabatic region of the neutrino mixing parameter space. The effect of the 
noise is large in the adiabatic region only and becomes increasingly negligible 
as the neutrino flavor transition becomes nonadiabatic. A 1% fluctuation will 
decrease Y e from 0.5 to ~ 0.495 for points in the MSW parameter space on 
the boundary of the excluded region [5]. However, we consider this rms 
fluctuation value to be actually too large for the relatively quiesent TPB > 3 
s epoch. For an order of magnitude smaller fluctuation which may be more 
reasonable, the effect of the noise is totally ignorable. In other words, noise 
at the 0.1% level will not alter the region of the neutrino mixing parameter 
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space excluded in Ref. [5] by the MSW effect. 

Noise with an rms amplitude of 0.05% of the averaged local matter density 
can lead to significant neutrino flavor decoherence during the shock reheating 
epoch at TPB ~ 0.15 s. In turn, this would lead to a dimunition of the 
adiabatic MSW-induced increase in the supernova explosion energy from ~ 
(30-60)% to ~ 20%. The physics of supernova explosion process is quite 
complicated and depends on many factors. In our opinion, if the heating 
of the stalled shock by an MSW transition is determined to be a necessary 
component for a successful explosion, then a more detailed analysis of the 
effects of random matter density fluctuations on the MSW transition would 
be warranted. 
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Table I. Ratios of the heating rates in the presence of MSW transitions 
with no fluctuations ((3 = 0), 0.05% fluctuations ((3 — 5 x 10 -4 ), and 0.1% 
fluctuations ((3 = 10 -3 ) to the heating rate in the absence of flavor transitions. 
The neutrino temperatures are T Ue = T Pe = 5 MeV and T Ur = 7 MeV. 
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(3 = 5 x 10" 4 


(3 = 10 
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io- 7 


1.31 
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1.19 


8 x 10 3 
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1.19 
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Figure Captions 



Figure 1. Survival probability P[y T — > v T ) is given as a function of the 
correlation length, r c for the numerical simulation by the method of Ref. [11] 
and by the approximation in Eq. (9). The parameters chosen were Sm 2 = 10 2 
eV 2 , sin 2 29 = 10" 3 , E = 33 MeV, and r c ° = L rcs /n = (5m 2 sin 29/ 'AE)- 1 = 
832.6 cm for (a). The values 5m 2 = 10 2 eV 2 , sin 2 29 = 10" 5 , E = 33 MeV, 
and r c ° = 8326 cm were chosen for (b). The rms value of the fluctuations 
was 1% of the local density in both (a) and (b). The matter density profile 
at TPB m 3 s was used [5] for both (a) and (b). The horizontal line in (b) is 
the survival probability in the absence of fluctuations. 

Figure 2. Survival probability P[y T — > v T ) is given as a function of energy 
for the numerical simulation by the method of Ref. [11] and by the approxi- 
mation in Eq. (9). The rms value of the fluctuations is taken to be 1% in (a) 
and 0.5% in (b). In both (a) and (b) the survival probability in the absence 
of the fluctuations is shown by the dashed line. The parameters chosen were 
5m 2 = 10 2 eV 2 , sin 2 29 = 10" 5 , r c = 8326 cm. The matter density profile at 
TPB w 3 s was used [5]. 

Figure 3. Survival probability P[y T — > u T ) is given as a function of energy 
for the numerical simulation by the method of Ref. [11] and by the approx- 
imation in Eq. (9). The rms value of the fluctuations is 0.05% in (a) and 
0.1% in (b). In both (a) and (b) the survival probability in the absence of the 
fluctuations is essentially zero. The parameters chosen were Am 2 = 10 3 eV 2 , 
sin 2 26* = 10~ 6 , r c = 2792.4 cm. The matter density profile at TPB w 0.15 s 
was used [4] and the density profile was cut off at a radius corresponding to 
the position of the stalled shock. 
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